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Abstract 

We present a short account of the present experimental situation of stripes in cuprates followed by a review of our 
present understanding of their ground state and excited state properties. Collective modes, the dynamical structure fac- 
tor, and the optical conductivity of stripes are computed using the time-dependent Gutzwiller approximation applied 
to realistic one band and three band Hubbard models, and are found to be in excellent agreement with experiment. On 
the other hand, experiments like angle-resolved photoemission and scanning tunneling microscopy show the coexis- 
tence of stripes at high energies with Fermi liquid quasiparticles at low energies. We show that a phenomenological 
model going beyond mean-field can reconcile this dynamic dichotomy. 

Keywords: high-temperature superconductors, stripes, electronic structure, transport properties, charge excitations, 
magnetic excitations 



1. Introduction 

Since the discovery of high-temperature supercon- 
ductors (HTSC) by Bednorz and Miiller [iT] numerous 
experiments have evidenced the existence of electronic 
inhomogeneities in these compounds (cf. e.g. Ref. iQl)- 
While early on these inhomogeneities where believed to 
be predominantly chemically driven, e.g. due to mate- 
rial imperfections like disorder induced by the dopant 
ions, it was subsequently realized that the strongly cor- 
related character of the cuprate superconductors and 
thus the electronic subsystem itself can favor the forma- 
tion of inhomogeneities on the nano-scale. According to 
the analyis by the Rome group jstl the inherent reduction 
of the quasiparticle kinetic energy together with a short 
range attractive (e.g. electron-lattice) interaction can in- 
duce an instability towards phase separation since com- 
peting nesting instabilities are suppressed due to a resid- 
ual repulsion of the quasiparticles at large momenta. On 
the other hand, the long-range repulsive Coulomb inter- 
action will spoil the associated zero-momentum insta- 
bility in the charge sector and instead shift the wave- 
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vector of the ordering transition to finite values. This so- 
called frustrated phase separation mechanism ^ thus 
results in an incommensurate charge ordering (CO). 



An alternative approach is based on Hartree-Fock 
(HF) investigations of Hubbard (and tJ)-type hamilto- 
nians [Hi 0, 13, Htl which have revealed solutions with a 
combined charge- and spin-density wave. In contrast 
to the frustrated phase separation mechanism the charge 
order in this case is driven by a spin-density wave in- 
stability which via the coupling between longitudinal 
spin and charge degrees of freedom results in a con- 
comitant charge density wave. These so-called stripe 
solutions have been confirmed later on by more sophis- 
ticated numerical methods. Within a density matrix 
renormalization group (DMRG) approach White and 
Scalapino 1 in have found stable domain wall so- 

lutions in the physically relevant doping regime of the 
tJ-model. The stripe stability in the tJ-model has also 
been investigated using exact diagonalization 11121 Il3n 
as well as quantum and variational Monte Carlo tech- 
niques HI 15i1 . Besides the HF approach a variety of 
methods has also been applied to Hubbard-type mod- 
els in order to investigate the possibility of charge and 
spin order. Fleck and collaborators [16] have found sta- 
ble stripes using dynamical mean-field theory (DMFT) 



Preprint submitted to Elsevier 



February 9, 2012 



and a cluster perturbation approach has been applied in 
Ref. 12^. From an ab-initio perspective stripe order in 
lanthanum cuprates has also been shown to be stable in 
LDA+U calculations lEHl . 

Here we present a short update of the current experi- 
mental situation (Sec. |2|i followed by an overview over 
static and dynamical properties of stripes ll22i l23l |25 , 
ll 111 Izlllllli, 111 H obtained within the time- 
dependent Gu tzwiller approximation (TDGA) il33i 
and phenomenological models. In Sec. 
Owe first derive the parameter set for the extended Hub- 
bard model which in the following sections allows us to 
make quantitative comparison with experimental data. 
We then show in Sect. |4]that the Gutzwiller approxi- 
mation (GA) of this model in fact leads to stable stripe 
ground states which allows for the doping dependent 
evaluation of the corresponding incommensurate spin 
response. In Sec. |5]we discuss the electronic structure 
and various transport properties of stripes. Sec. |6]fo- 
cusses on the momentum and frequency dependent spin 
excitations where we also address recent resonant in- 
elastic x-ray scattering experiments. Sec. |7]presents re- 
sults for the charge excitations from stripe ground states, 
focussing on the optical conductivity. Finally, we show 
in Sec. [8] how our Gutzwiller variational theory can be 
extended towards a dynamical description of stripes be- 
fore we summarize our discussion in Sec. |9l 

In the following section we start with an update of the 
experimental situation. 



2. Brief Experimental Review 

First evidence for stripe order in HTSC's came 
from elastic neutron scattering experiments by Tran- 
quada and collaborators 131. 32L 13311 . They observed 
a splitting of both spin and charge order peaks in 
Lai.48Ndo.4Sro.i2Cu04 (LNSCO) which resembled sim- 
ilar data in the nickelates where both incommensurate 
antiferromagnetic (AF) order lf34l [36ll and the ordering 
of charges ll35i[36ll has been detected by neutron scatter- 
ing and electron diffraction, respectively. In Ref. ll36ll 
it was shown that the magnetic ordering displays itself 
as an occurence of first and third harmonic Bragg peaks 
whereas the charge ordering is associated with second 
harmonic peaks. Meanwhile analogous static stripe or- 
der has also been detected in Eu -codoped lanthanum 
cuprates (LECO) |38] and La2-,Ba,Cu04 (LBCO) |i37» 
dOII which all display a LTT lattice distortion. That the 
latter is not necessarily a condition precedent to the for- 
mation of stripes has recently been demonstrated by ap- 
plying pressure to a Lai 875Bao.i25Cu04 sample. This 



restores the structure to fourfold lattice symmetry keep- 
ing symmetry broken electronic stripe states 1 39]. More- 
over, charge and spin stripe order can also be induced in 
the LTO phase of im puri ty (Cu, Zn, Fe, Ga) substituted 
lanthanum cuprates [40] without the need of applying 
pressure. Whereas charge order in the neutron and hard 
x-ray scattering 114 ill experiments mentioned above has 
only been detected indirectly via the associated lattice 
modulation, more recent soft resonant x-ray scattering 
studies on LBCO IH and LECO H have directly 
revealed the spatial modulation of charge in these sys- 
tems. 

The incommensurability 6 (defined as the shift of 
the magnetic Bragg peaks from the antiferromagnetic 
(AF) wave-vector Qaf in reciprocal lattice units) in 
samples with long range stripe order follows the so- 
called ' Yamada-Plot' 14411 and depends linearly on dop- 
ing 6 = «/, up to ni, X 1/8. Here «/, is the number 
of added holes per planar Cu with respect to the par- 
ent insulating compound. This behavior is compatible 
with stripe like modulations of charge and spin having 
a linear concentration of added holes v = 1 /2 (so-called 
half-filled stripes) as we will discuss in more detail in 
Sec.m 

Evidence for some kind of stripe order in other HTSC 
materials is based on the doping dependence of the 
low energy spin response. In non-codoped lanthanum 
cuprates (LCO) the incommensurability 6 (defined in 
terms of the shift of the low energy magnetic scatter- 
ing from Qaf) also follows the 'Yamada-Plot' of the 
samples with long-range order. This suggests that the 
static stripe order is replaced by some kind of 'fluctuat- 
ing order' in the Nd- (or Ba-, Eu-) free samples. Above 
«/, ~ 1/8, the incommensurability stays essentially con- 
stant but the intensity of the low energy spin fluctua- 
tions decreases and vanishes at the same concentration 
where superconductivity disappears in the overdoped 
regime [45]. This behavior is mirrored at low doping in 
that superconductivity, in LSCO materials, disappears 
upon decreasing doping exactly at the point in which 
incommensurate scattering parallel to the CuO bond 
disappears. |49| This strong correlation between super- 
conductivity and low energy incommensurate scattering 
parallel to the CuO bond suggesting an intimate relation 
between both phenomena and provides a strong motiva- 
tion to understand stripe physics in cuprates. 

Once superconductivity disappears by decreasing 
doping below «/, ^ 0.055 in lanthanum cuprates in- 
commensurate scatterin g do es not disappear entirely but 
rotates by 45° S 113 lH, 51] to the diagonal direc- 
tion. Additionally the orthorhombic lattice distortion 
allows one to conclude that the elastic diagonal mag- 
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netic scattering is one-dimensional with the associated 
modulation along the orthorhombic b*-a\is, supporting 
again the picture of stripe formation. When S is mea- 
sured in units of reciprocal tetragonal lattice units in 
both the vertical and diagonal phase it turns out that 
the magnitude of the incommensurability numerically 
coincides across the rotation leading to a linear rela- 
tion 6 - rij,. Upon approaching the border of the AF 
phase at «/, = 0.02 the incommensurability approaches 
e - «/, [51] where e is measured in units of reciprocal 
orthorhombic lattice units, thus e = The crossover 
from diagonal to vertical incommensurate spin response 
is so far only observed in lanthanum cuprates whereas 
in all other HTSC compounds the scattering is always 
along the Cu-O bond direction. Despite this peculiarity 
of LCO other experiments point to universal behavior. 

A linear relationship between doping and incom- 
mensurability is also observed in YBCOe+A although 
the curve falls below the one of lanthanum cuprates 
I 52II53II54II . Apart from details, this suggests that stripe 
physics is a universal phenomenon in cuprates. 

Upon increasing frequency of spin excitations the 
dispersion develops a spectrum ('hour glass') which is 
similar in many high-T^ cuprates again suggesting uni- 
versaUty. At a certain frequency £5 the incommensu- 
rate branches merge [55] at the AF wave-vector and for 
even higher energies the excitations become incommen- 
surate again. The 'hour glass' sgectrum has been de- 
tected in the vertical fs^, 58, 5^ 63l and diagonal phase 
1 62L 164. ] of lanthanum cuprates, the YBCO compounds 



codoped lanthanum cuprates, YBCO and bimuthates. 
Experimental data of local probes like NQR (sE, 1^ 



Ml, NMR [88], EXAFS [182], and X-ray microdiffrac- 




and Bi2Sr2CaCu208+v 1170117 lH 



Magnetic fluctuations which dispersion follows the 
'hour-glass' shape naturally arise from stripe correla- 
tions in the ground state Izi S IM S S H S [tI, 



M m m l83l |29t[ (cf. Sec. ^. BasicaUy a Gold- 
stone mode arises from each of the incommensurate 
wave-vectors and disperses in a cone-shaped structure 
to higher energies. However, due to disorder and pecu- 
liar magnetic couplings across the stripes (as supported 
by recent inelastic neutron scattering (INS) experiments 
on La2-.iSrvCo04 [^8?]) the weight of the outwards dis- 
persing contribution is strongly suppressed so that the 
excitations which dominate the spectral weight disperse 
towards the AF wave-vector Qaf- The corresponding 
excitation energy Eg at Qaf is essentially determined by 
the spin coupling across the antiphase domains. Since 
at high energies the excitations should again resemble 
those of the (undoped) AF one finds an outwards dis- 
persing intensity also above Eg. 

While apparently the spectrum of spin excitations 
is compatible with stripe correlations the problem is 
the detection of an associated charge order in the non- 



tion 1118311 are suggestive of the formation of electronic 
inhomogeneities, possibly induced by magnetic field 
isol] . but cannot provide evidence for long-range order. 
In this context it is remarkable that due to refinements 



in the experimental technique Haase et al. 18811 where 
able to demonstrate a correlation of charge and density 
variations on short length scales. 

One way to reconcile the 'stripe-like' magnetic exci- 
tation spectrum with the absence of long-range charge 
order is the notion of a charge nematic [90] which 
(in the charge channel) breaks rotational but not trans- 
lational invariance. This concept is extremely use- 
ful in detwinned YBCO where both neutron scattering 
168, 169I1 and thermoelectric transport (9?] data show 
pronounced in-plane anisotropics but may also apply to 
the spin-glass phase of lanthanum cuprates where, as 
mentioned above, the incommensurate spin scattering 
is one-dimensional. 

Due to the absence of clear signatures of charge or- 
dering in other than co-doped lanthanum cuprates alter- 
native theories have been proposed in order to account 
for the incommensurate spin scattering and the associ- 
ated 'hourglass' magnetic spectrum ^92, 93, 94, 95, 9^ 
97, 98, 99. ,100]. Nonetheless there are a number of 
fingerprints in spectroscopic probes which support the 
presence of incommensurate charge scattering. For ex- 
ample, the anomalous lineshape and temperature depen- 
dence of the bond-str etchin g phonons in various cuprate 
materials (cf. Ref. [^ and references therein) can 
be attributed to the influence of stripes. For example, it 
has been proposed that the softening of these phonons at 
certain momenta is due to the influence o f cha rge stripe 
fluctuations on the phonon self-energy il02ll In addi- 
tion one may also have a Kohn-typ e ano maly due to the 
nesting along the half-filled stripe Furtheron, the 

scattering of charge carriers by incommensurate charge 
fluctuations shou l d ind uce a reconstruction of the Fermi 
surface 1 104 1 10511 ■ Angle-resolved photoemission 
experiments (ARPES) on LSCO [ 106, 107] have indeed 
shown the appearance of straight Fermi surface (FS) 
segments around the M-points of the Brillouin zone as 
expected for a striped ground state. However, one pecu- 
liar feature of these experiments concerns the fact that 
the stripe-like FS was obtained from the momentum dis- 
tribution nii, by integrating the spectral function over a 
broad energy window (~ 300 meV). On the other hand, 
upon following the momentum dependence of the low- 
energy part of the energy distribution curves a large FS 
was found corresponding to the LDA band-structure and 
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fulfilling Luttinger's theorem. We will address this di- 
chotomy of stripe-like spectral features at large energies 
and a 'protected' FS at low energies in Sec. [8] 

In LECO, where static charge order has been de- 
tected with RIXS [43'], the associated reconstruction 
of the FS has indeed been resolved by ARPES. More- 
over, a recent series of magnetotransport experiments 
I lld.llllllll2 , 113 ] have revealed quantum oscillations 
suggestive of the formation of electron pocktes due to 
a FS rec onstr uction. It has been shown by Millis and 
Norman fTTT] that a charge density wave ground state 
with intermediate values of the stripe order parameter 
is in qualitative agreement with the quantum oscilla- 
tion data. A recent comparative study of thermoelec- 
tric properties in YBCO and LECO 1 1 15] also provides 
strong evidence that the Fermi surface reconstruction in 
YBCO is associated with stripe correlations. In the con- 
text of transport experiments one should also mention 
the Nernst effect which in hole doped HTSC 's above 
usually yields a large positive signal [116] and which 
has been considered as indicative of fluctuating super- 
conductivity. However, in addition also a FS reconstruc- 
tion contributes to the Nernst signal and both contribu- 
tions have been resolved in LNCO and LECO I iTtIi . 
These authors have attributed the high temperature sig- 
nal to originate from changes in the FS, an interpreta- 
tion which has recently been challenged in Ref. [iTsll . 
A theoretical analysis of the Ner nst effect in terms of 
stripe order can be found in Refs. ifTHITIoll . 

Scanning tunneling microscopy (STM) provides ad- 
ditional evidence (limited obviously to the surface) for 
charge order in cuprate superconductors. Correspond- 
ing measurements performed on bismuthate and oxy- 
chloride superconductors have revealed a complex mod- 
ulation of the local density o f states (LDOS) both i n the 



supe rconducting (S C) state 1 121 



122, 123, 124, 125 



126 1 and above Tc 1 126l 127 , 128 1. In both cases one 
observes peaks in the Fourier transform of the real space 
LDOS at wave-vectors Q - 27r/(4ao)---27r/(5ao) sugges- 
tive of checkerboard or stripe charge order However, it 
is important to distinguish between the case where these 
peaks are non-dispersive in energy (and thus signature 
of 'real' charge order) or the situation where they foUow 
a bias-dependent dispersion due to quasiparticle (QP) 
interference. In the latter case the spatial LDOS varia- 
tions can be understood from the so-called octet model 
lfl2ll[l29^ which attributes the modulations to the elas- 
tic scattering between the high density regions of the 
Bogoljubov 'bananas' in the superc onducting state. Re- 
cent STM investigations lfl30[[l3?l may resolve this ap- 
parent conflict since they suggest that both, dispersive 
and non-dispersive scattering originates from diff'erent 



regions in momentum and energy space. The states 
in the nodal region which are well defined in k-space 
and undergo a transition to a li-wave SC state below 
Tc are then responsible for the low energy QP interfer- 
ence structure of the LDOS, whereas the ill-defined k- 
space 'quasiparticle' states in the antinodal regions are 
responsible for the non-dispersive CO above some en- 
ergy scale (jL>Q. As mentioned before we will address the 
issue of high energy electronic order in Sec. |8] 



3. Model and parameters 

Moment formation in a spin-density wave solution in- 
volves the large Hubbard U scale and modifies the elec- 
tronic structure of the system over a large energy range. 
Therefore, it is not a priori obvious that the reduction of 
a multiband model to a more simplified model hamilto- 
nian will retain the relevant features in the spectra, es- 
pecially when one is concerned with high frequency ex- 
citations. A suitable method in order to deal with this 
problem has to fulfill essentially the following two con- 
flicting criteria: faj it should be sufficiently accurate in 
order to account for the Hubbard-type correlations of 
the model and (b)\\. should allow for the investigations 
of large systems in order to provide a sufficient resolu- 
tion for the correlation functions in momentum space. 

Subsequent to the first stripe calculations JHQ^S^ 
investigations of collective excitations have been per- 
formed within the Hartree-Fock (HE) approximation of 
the one-ban d Hubbard model supple mented with RPA 
fluctuations |[i33,[iM[j39',T40', W]. However, in case 
of the HE approximation] 142] it turns out that the stripe 
solutions relevant for cuprates are only favored for un- 
realistic small values of U It ^ 3... 5 whereas a ratio of 
f//f a; 8 is required to reproduce the magnon spectrum 
ofLCO (cf. below). 

In order to enable a quantitative comparison with ex- 
periment, the results we will discuss in the following 
sections will b e mos tly based on the Gutzwiller approx- 
imation (GA) 1 132! supplemented with Gaussian fluc- 
tuations II 33, . 134,.. 135... 136.1 (so-called time-dependent 
GA [TDGA]). 

We consider the one-band Hubbard model 



H 



z 



(1) 



where c,> (cj^) destroys (creates) an electron with spin 

o" at site /, and n,-o- - c'. c^j. U is the on-site Hub- 
bard repulsion and t/j denotes the hopping parameter 
between sites / and j. We restrict to hopping between 
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nearest (~ f) and next-nearest (~ f')neighbors. Forhigh- 
Tc cuprates the value of U can be estimated fro m th e 
magnon dispersion of the undoped system [ 144l 1 14511 . 
In fact, within spin-wave theory (SWT) applied to the 
Heisenberg model (corresponding to f//f — > oo) the 
magnon excitations are given by 



-^1 - [cos + cos ^y]^/4 



(2) 



and thus the dispersion vanishes along the magnetic 
Brillouin zone. Including quantum fluctuations beyond 
SWT [143i] leads to excitations with lower energy at 
(tt, 0) than at (7t/2,7t/2), contrary to wha t is o bserved 
in undoped cuprate superconductors 111441 114511 . It has 
been argued that in cuprates corrections to the Heisen- 
berg model a rising as higher orders in a f/I/ expansion 
are relevant lll46l ll47l Il48i [Hi, [lH. The most im- 
portant of such corrections is a term which cyclically 
exchanges four spins on a plaquette. A sizable value 
for this term has been revealed by analyzing phonon- 
assisted multimagnon infrared absorp tion! 1491 and the 
dispersion relation measured with INS| 144 , 145ll shown 
in Fig. [1] In particular, the dispersion in the (n, 0) and 
{7t/2, n/2) is mainly due to this term. In case of the Hub- 
bard model the magnetic zone edge dispersion is due 
to the finiteness of U/t. For moderate values of U/t 
this leads to dispersion which decreases from {n, 0) to 
{ji 12,71 12) along the magnetic zone boundary whereas 
for large Ujt the behavior of the spin- 1/2 Heisenberg 
model is recovered l,1501 . 

Fig. [T^ shows the J//f dependence of the magnetic 
excitations for t he ha lf-filled AF as compared to INS 
data from Ref. 1 14411 . The dispersion along the mag- 
netic zone boundary can be accurately fitted by a value 
of Ujt - 8 where the absolute energy scale is fixed by 
the nearest neighbor hopping t - 36QmeV. Panel b of 
Fig. [1] displays the dependence of the spin excitations 
on the next-nearest neighbor hopping t' /t. For small to 
moderate values of f'/f the zone boundary dispersion is 
only slightly dependent on this parameter, however, sig- 
nificant softening occurs for f'/f = -0.5. The increas- 
ing frustration between nearest and next-nearest neigh- 
bor magnetic interactions drives the system towards an 
ID-AF instability the precursor of which is seen as a 
magnetic mode softening at q = (n, 0). 

Further information on the parameter f'/f can be ob- 
tained from the optical conductivity which in the inset 
to Fig. [TJ) is shown for the same parameters as used in 
the main panel. The lowest energy excitation Q„„„ for 
the half-filled system is shifted up with increasing \t' /t\ 
and the ex perimental value of 0,„,„ x 2.leV for LCO 
1 15ll I152II is accurately reproduced for the parameter 
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Figure 1 : Energy and wave vector dependence of magnetic excitations 
in tlie lialf- filled system as obtained from the TDGA for (a) f /t = 
-0.2 and varying U/t; (b) U/t = 8 and varying t' /t. Squares and 
triangles correspond to data points from INS experiments on LaCu04 
by Coldea et al. 1144], The inset to (b) shows the optical conductivity 
for the same parameters as used in the main panel (b). 



set U/t = 8, t'/t = -0.2, t = 360meV. In the follow- 
ing section we will show that a value of f'/f = ~0.2 is 
also appropriate to account for the doping dependent in- 
commensurability of lanthanum cuprates 14411 based on 
a striped ground state. 

Note that some few results discussed in the present 
review (Sees. 14.11 1711 are also obtained within the three- 
band model which for copper li.v^-v^ ™d oxygen p^ y 
orbitals includes the associated hopping processes, or- 
bital energies, and intra- and interorbital Coulomb re- 
pulsions. We have taken these parameters either from 
first principle computati ons [167. .168.1 or from more 
empirical considerations il56ll . 
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4. Stripe stability and static properties 



4.1. Charge and spin structure 

Stripes are characterized by one-dimensional an- 
tiphase boundaries for the AF order which also host 
the doped holes. These textures can be understood as 
the strong coupling remna nt of an instability in the spin 
channel. It can be shown 1 154'] that for dispersion rela- 
tions relevant for cuprates the orientation is either along 
the copper-oxygen (vertical) or the diagonal direction. 
The domain walls can either reside on lattice sites or 
bonds and these two prototypical stripe solutions for 
both vertical and diagonal orientations are sketched in 
Fig. El 

By symmetry, for site-centered (SC) stripes the mag- 
netization vanishes on the domain wall whereas for 
bond-centered (BC) stripes these are build up from fer- 
romagnetically aligned nearest-neighbor spins. For di- 
agonal BC stripes this leads to the situation that con- 
figurations with d - even acquire a macroscopic mag- 
netization in contrast to vertical BC stripes where the 
ferromagnetic bonds alternate along the domain wall. 
Diagonal BC stripes can be viewed as a realization of 
the staircase structures proposed by Granath I.155J with 
step length 1-2. 

Fig. [3] shows the charge and spin density for d - 
4 BC vertical stripes (doping «/, - 1/8 ) obtained 
withi n the 3-band model and parameters taken from 
Ref. lll56ll . A similar charge distribution was pre- 
dicted in Ref. lESt] and found to be in excellent agree- 
ment with a charge sensitive probe by Abbamonte and 
collaborators I'^S]. However, these authors could not de- 
termine if the stripes where Cu centered or O centered. 
More recently Davis and coUaborators lI 1 5711 have im- 
aged glassy stripes which indeed are centered on O as 
predictedl 25| . Taken together these experiments give us 
amplitude and phase information of stripes in excellent 
agreement with Ref. f23\ thus showing that it is possi- 
ble to obtain and even predict realistic information on 
the intermediate scale physics of cuprates. 

The Fourier transform of the charge (oc pg) and spin 
distribution (oc niQ) determines Bragg peak weights 
(oc p^, OTg ) in s cattering experiments (for the definitions 
see Ref. fl58^) as shown in Fig. |4l Disregarding the 
difference in cross section for different processes (mag- 
netic neutron scattering vs. X-ray or nuclear neutron 
scattering) we see that Bragg weights for the charge are 
practically 3 orders of magnitude smaller than those for 
the Cu spins. This is due to very soft charge distribu- 
tion shown in Fig. [3] with respect to the spin distribution 
and explains why it has been so hard to detect charge 
ordering. 



M i t 

i t \ 

t i t 

i t \ 

t i t 



b) I I 

i t 

t i 

i t 

t i 



t 



t t 
t t 
t t 



i t 
t i 
i t 
t i 
i t 

i t 
t i 
i t 
t i 
i t 





dtet 









t i t 




t i 




1 t i 


t 


i t 


t 


\t i t 


i t 




i t 




t i 1 1- 




* t t i 


t \ t 


1 1 1 i 


t i 


it t 


1 t i 


t i 1 1 


i t 


1 1 1 


t i t 


i t i t 


t i 


i t i 


t t \ 


t i t i 


t t 







dtet 





Figure 2: Spin structure of vertical (a,b) and diagonal stripe textures. 
The doped holes mainly reside on the (shaded) domain walls. In 
both cases the prototypical structures comprise site-centered (a,c) and 
bond-centered structures (b,d). The distance between charge stripes, 
measured along the copper-oxygen bond direction is denoted as d in 
units of the tetragonal lattice constant a,,,,. 
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Figure 3: Charge and spin density for d = 4 BC vertical stripes. Dop- 
ing rill = 1/8. Open circles (p orbitals) represent Cu (0) sites. The 
numbers and the size of symbols represent the excess charge whereas 
the spin density is proportional to the length of arrows. 
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Figure 4: Bragg weights m^, at reciprocal lattice vectors (indicated 
by the dots) for «;, = 1/8 and d = 4 stripes. The polygons around each 
point are the reduced magnetic Brillouin zones. We use reciprocal lat- 
tice units of the fundamental Brillouin zone. O,- (Oy) are the oxygens 
in the bonds perpendicular (parallel) to the stripe. 



4.2. Stripe stability and incommensurability 

The task is now for a given parameter set and dop- 
ing to variationally determine the lowest energy state 
among the various stripe states and to compare with 
other possible solutions. This has been accomplished 
in Refs. |27, isijsoll for the one-band hamiltonian Eq. 
[T]and in Ref. |25] also for the three-band model. Fig. |5] 
displays the binding energy per hole 

EAFiNh = 0) - E,M,) 



Nh 



(3) 



for the energetically most competing structures, namely 
stripes, spirals and spin polarons. Here Eaf(Ni, = 0) 
denotes the energy of the undoped AF and E,ext{Nh) 
is the energy of a given texture obtained for doping 
the system with Nh holes. We can also define the fill- 
ing factor V of stripes which corresponds to the con- 
centration of holes per unit cell along the stripe. This 



can be related to the hole concentration n^ and the dis- 
tance (in lattice units) between charge stripes d as 
V - ni, ■ d. From the inset to Fig. [5^ it can be seen that 
for U/t - ^ and t' /t - -0.2 low doping vertical stripes 
have Vopi ~ 0.55 whereas the minimum of diagonal SC 
(BC) domain walls is at Vopr = 1(0.75). We will see be- 
low that Vopt determines the slope of the Yamada plot at 
low doping. 

We find that for the one-band model, vertical BC 
and SC stripes are also practically dege nerate in energy 
in agreement with Refs. |@, [3, 163 1. For definite- 
ness we mainly restrict ourself to the BC case because 
these textures constitute the more stable configuration 
at «/, = 1/8 in the more accurate three-band model lEsIl 
and in first principle computations ^ ill . However, one 
should keep in mind that all energetic considerations be- 
low hold equally for SC vertical stripes. 

From the inset to Fig. |5^ it turns out that diagonal 
BC stripes are practically (accidentally) degenerate in 
energy with vertical stripes. The energy of the DSC tex- 
ture is ^ 0.02f per hole above. These small energy dif- 
ferences should not be significant given the simplicity of 
the model. A precise determination of the relative sta- 
bility of the different phases would require at least mul- 
tiorbital effects, inclusion of both long-range Coulomb 
interactions and coupling of the holes to the tilts of the 
Cu04 octahedra. The latter have been shown to play a 
major role in the stabilization of vertical vs. diagonal 
stripes flS?]. 

In Fig. |5]the dot-dashed line corresponds to the en- 
velope of the stripe binding energy curves for different 
periodicites d. The solid (dashed) lines show e/, for di- 
agonal (vertical) spiral solutions which apparently are 
unstable towards phase separation at low doping. The 
horizontal line is the energy of the phase separated state 
obtained from the Maxwell construction. The respec- 
tive stability of stripes and uniform spirals is determined 
by the ratio between next- nearest ne ighbor and nearest 
neighbor hopping t'/t jsoi [l6ll [l62ll . For t'/t = -0.2 
stripes are stable over the whole doping range up to 
ni, ^ 0.28 where their energy becomes degenerate with 
that of vertical spirals. In contrast, it is found that with 
increasing \t' /t\ stripes with significant spin canting start 
to dominate the phase diagram and for even larger \ t' /t\ 
also checkerboard textures may become the ground state 
128. 30ll . Also the optimum filling Vopi changes with t'/t 
as expected according to the results of Ref. frf]. The 
influence of a next-nearest neighbor hopping term on 
the stripe stability has been previously investigated by 
various methods. From the DMRG approach applied to 
the t-t'-J model I IJJ it turned out that a negative t' /tcan 
suppress both the formation of stripes and pairing cor- 
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Figure 5: (a) Binding energy per hole e;, as a function of doping for 
stripes, uniform spirals and spin polaron. The inset shows e/, for d = 
10 vertical (dashed-dotted), bond centered diagonal (dashed), and site- 
centered diagonal (solid) stripes as a function of the fining factor v. (b) 
Incommensurability 1 / (2d) of vertical stripes and uniform spirals as 
a function of doping. Diamonds are data for LSCO from Ref (4^; 
full dots are data for Fe-codoped LSCO from Ref. 1 160]. Results 
are obtained from the one-band model with parameters: U/t = 8 and 
t'/t = -0.2. 



relations. A weakening of stripe tendencies for t'/t < 
in the same model was also found with exact diago- 
naUzation lll3ll and in the Hubbard model with DMFT 
lUS], and the HF approximation. iHiiiliil AH these 
calculations suggest that static stripes are destabilized 
when the ratio f'/f becomes negative. This may indicate 
a more dynamical character of stripes in some systems 
like Tl and Hg based compounds. 

With regard to the experimental evidence of stripe 
correlations in lanthanum cuprates it turns out that a 
ratio of f'/f = -0.2 can account for the doping de- 
pendent incommensurability 6 observed in these com- 
pounds [44]. For vertical stripes this quantity is related 
to the modulation of the AF order and thus given by 



(4) 



in units of 2nla,^, where a,e, is the tetragonal lattice con- 
stant. The charged core of the stripe has a characteristic 
width f so that when the charge periodicity d is larger 



than ^ there are negligible interstripe interactions and 
doping proceeds by increasing the number of stripes. 
In this regime the number of stripes for a fixed number 
of holes is optimized when v = Vopt- Thus Vopt deter- 
mines the slope of the incommensurability as a function 
of doping, according to Eq. (|4]i and Vopi - 0.55 implies 
6 - 0.91n/, in good agreement with the Yamada plot. 
The regime d > ^ is, characterized by the fact that e/, at 
the minimum is independent of doping or d. From Fig.|5] 
we find ^ - A- lattice units which coincides with a direct 
examination of the charge profile. For d < ^ stripes 
completely cover the plane and increasing further the 
number of stripes becomes energetically costly. There- 
fore doping proceeds by increasing the charge of stripes 
and 6 becomes constant. The crossover in doping be- 
tween the two regimes occurs when n/, = Vop,/d ~ 1/8 
in good agreement with the Yamada plot in LCO. 

As can be seen from Fig. |5}' the incommensurability 
for stripes shows a staircase structure where the steps 
occur at the crossing of the corresponding energy curves 
(cf. Fig. |5^). For small doping (d > since interstripe 
interactions are negligible, one can produce a practically 
continuous curve by considering combinations of solu- 
tions with periodicity d and d + I. As explained above 
this is not any more convenient for d < ^ and the ad- 
ditional holes are doped into existing stripes producing 
the plateau in the incommensurability. 

Contrary to the stripes, one observes a continuous in- 
crease of (5(«/,) in case of uniform spirals . Neverthe- 
less, both phases show similar evolution of incommen- 
surability with doping and indeed fairly similar values 
of 6(ni,) up to «/, ~ 1/8. For larger doping this agree- 
ment is lost due to the tendency of stripes to make wide 
plateaus. 

We have explained the Yamada plot in terms of well 
formed stripe textures. It is remarkable that the incom- 
mensurabilities, both for stripes and spirals, are not far 
from what one would obtain from an analysis of the 
momentu m dep endent susceptibility of the Fermi liq- 
uid phase 11531]. In fact, if one adiabatically decreases 
the interaction both uniform spirals an d str ipes origi- 
nate from the same magnetic instabilitv M 15411 . One can 
call this a "weak coupling" instability. On the other 
hand within the time-dependent Gutzwiller approxima- 
tion, the Hubbard interaction U is strongly screened by 
vertex corrections. Therefore a quite large value of U 
is needed in order to make the Fermi liquid unstable 
and the instability is better characterized as "interme- 
diate coupling". 

The time dependent Gutzwiller analysis shows that it 
is the plateau in the spin susceptibility close to Qaf - 
{n, n) (cf. Fig. 1 in Ref. lll54ll ') which drives the system 
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unstable towards spiral or stripe order for sufficiently 
large on-site repulsion Ujt. 

Figure also shows d ata fr om more recent neu- 
tron scattering experiments lll60ll on Fe-LSCO. In these 
overdoped samples an elastic incommensurate spin re- 
sponse was found close to the dominant nesting vectors 
as extracted from ARPES experiments. It was thus con- 
cluded that the induced incommensurate response sig- 
nals an inherent instability of the itinerant charge carri- 
ers, being different from the low doping stripes arising 
from localized Cu spins. However, from Fig. |5j) it turns 
out that the data are rather close to the incommensura- 
bility curve of stripes obtained for the LSCO parameter 
set. Remarkably this behavior was predicted in Ref . lEsIl 
much before the measurement became available. Since 
these stripe calculations are based on an itinerant ap- 
proach, there is no 'dichotomy' between low doping 'lo- 
calized' spin stripes and large doping itinerant ones, but 
both appear as different limits of the same model. 

Finally we briefly address the incommensurability of 
diagonal stripes. In tetragonal units the incommensura- 
bility is a factor of V2 larger than that for vertical stripes 
given in Eq. (HJi. The more stable BC diagonal textures 
have Vop, ~ 0.75 so that 6^iag = «/,/( V2v„^,) » n/,/1.06. 
Therefore the doping dependence of 6 for BC diago- 
nal and vertical stripes is almost identical in agreement 
with neutron scattering data [46,,47„.48, 491. One should 
keep in mind the effect of disorder which in the diago- 
nal phase reduces the correlation to the same order or 
even smaller than the periodicitv ll47ll whereas for ver- 
tical stripes the correlation length can reach 150fl,f, or 
around 20 times the magnetic stripe periodicity. This 
becomes especially important in context of spin excita- 
tions from stripes in the spin glass phase (cf. Sec. |6]l. 

5. Electronic structure and transport properties 

Fig. |6] shows the bandstructure for half-filled d - A 
bond-centered stripes which are oriented along the y- 
direction. The electronic states of the undoped AF or- 
dered regions (magnetization per spin m) give rise to 
lower and upper Hubbard bands (LHB/UHB) which are 
separated by A ~ f/ ■ w. Stripes induce the formation of 
additional bands in the gap and for partially filled stripes 
(i.e. < V < 1) one of them ('active band') crosses the 
Fermi energy. For the present case of 'half-filled' stripes 
the crossing occurs at k^. - +kIA for k^t - {a,et =1). 
Perpendicular to the domain walls the bands are rather 
flat with the residual small dispersion due to the overlap 
between adjacent stripes. 

The 'active band' implies a quasi one-dimensional 
Fermi surface at momentum ky - jt/4 (cf. Fig. |7^). 
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Figure 6: Band structure for d = 4 bond-centered stripes oriented 
along the y-direction. Tlie dispersion is sliown for a cut along the 
stripe and ky = 0. The electronically relevant active band (AB) crosses 
the Fermi energy (doping ni, = 0. 125) and is in the gap between lower 
(LHB) and upper (UHB) Hubbard band. The AF modulation in the 
direction of the stripes leads to a doubling of the peiiodicity so that 
< ky < jt/2. 



The corresponding weight in the full Brillouin zone is 
concentrated around the (tt, 0) point which results in 
the two parallel FS segments. There is also a cross- 
ing at ky = 37r/4 with less spectral weight and the 
associated segment is shifted more towards k^ = 0. 
In LSCO one expects a superposition from x- and y- 
oriented stripes which results in the FS shown in Fig. 
iTj). Such a prediction is in good agreement with ARPES 
data 1 106, 107] and even fine details as the small curva- 
ture of the segments due to interstripe hopping are re- 
produced. It sh ould be mentioned that the experimental 
data [10611 110711 strongly depend on the energy window 
which is chosen to extract the FS. We will address this 
point in Sec. [T] 

As discussed above, stripes are approximately half- 
filled for doping < 1/8 and for n/, > 1/8 additional 
holes are doped into the active band. Therefore the 
chemical potential is expected to be approximately con- 
stant for «/, < 1 /8 and decreases for «/, > 1 / 8 in qualita - 
tive agreement with the observed behaviorOlO 



109D 



The rate of change of /i with doping, being a high 
derivative of the energy, is very sensitive to finite size ef- 
fects and, moreover, few experimental points are avail- 
able in this doping range in order to allow for a precise 
comparison. A rough estimate indicates that the theo- 
retical rate of change of fi with doping for «/, > 1/8 is 
approx imately a factor of 2 larger than the experimen- 
tal one lll08 , 109 1. This may be attributed to an under- 
estimation of the mass renormalization in mean-field. 
Another possibility which goes in the right direction is 
phase separation among the c/ = 4 stripe solution and 
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Figure 7: Fermi surface for d = 4 bond centered stripes at n/, = 0.125. 
(a) y-axis oriented stripes; (b) superposition of x- and y-axis oriented 
stripes. 



the paramagnetic overdoped Fermi liquid which is also 
suggested by the doping evolution of magnetic excita- 
tions (cf. Sec. |6]l. 

We finally note that also anomalies in the Hall and 
Nernst coeflicient can be attributed to the FS re construc- 
tion due to striped ground states lf25lll 19i[l20ll . 



6. Spin excitations 
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Figure 8: Constant frequency cuts of the magnetic excitation spec- 
trum at (J = 5,55, and l05meV for the d = 4,ni, = 0.125 BC stripe 
structure. The momentum space corresponds to the magnetic Bril- 
louin zone which is rotated by 45" with respect to the stripe direction 
(cf. Fig.|9). 



Fig. [8] reports constant frequency scans of the 
imaginary part of the transverse magnetic susceptibil- 
ity = -/»4 / e'"'{TS:,(t)S%m for y-axis ori- 
ented d - A BC stripes in the magnetic Brillouin zone 
(cf. Fig. ID. 
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Figure 9: Subdivision of tiie full Brillouin zone into reduced zones 
defined from the elementary cell of y-axis oiiented d = 4 BC stiipes. 
Arrows (a) ... (d) indicate the scan direction for the magnetic disper- 
sions shown in Figs. ll0l and ll2l The intensity plots of Fig.[8]are shown 
within the magnetic Brillouin zone indicated by the dashed square. 



The low frequency intensity is concentrated around 
the wave vectors - (7r(l ± 1/4, tt) corresponding 
to the spin periodicity of the underlying d - 4 stripe 
structure. In fact, it turns out that the weight of the 
elastic magnetic peaks is much larger at Qs than at the 
higher harmonics (cf. Fig. |4] where the corresponding 
weights obtained from the 3-band model calculations 
are shown). Note also that the magnetic weight vanishes 
identically on the 'nuclear' Bragg points 2« - n(^/2, 0) 
which only contribute in the charge channel. As a result 
the dispersion along the cut (q^, n) (cf. Fig. [TOl «/, = 
0.125) shows the expected behavior for spin-waves in a 
striped system S [sl [HB ■ Starting from the Gold- 
stone mode at Qs one observes two branches of spin 
waves where the large intensity one disperses towards 
Qaf = The other branch rapidly looses intensity 

which is analogous to what is found within linear spin- 
wave theory of the Heisenberg model 1731 ISll l82l 18311 
for small ratios of J_i_/J\\. Here J± denotes the (ferro- 
magnetic) exchange coupling across a stripe whereas J\\ 
denotes the coupling within the AF ordered legs. In ad- 
dition the intensity of the outwards dispersing branch is 
strongly reduced by disorder Iii84il . 

At £■ a; 55meV the spin excitations have finally 
shifted towards the AF wave-vector producing an in- 
tense feature at Qaf (middle panel of Fig. [8]l which 
is associated with a saddle-point in the magnetic dis- 
persion [Fig. [TO](n/, - 0.125)]. At the same time one 



observes the appearance of intensity in the stripe direc- 
tion starting from Qaf- From Fig. [TO](n/, - 0.125) it 
turns out that this feature is due to a slight softening of 
the magnetic excitations upon dispersing away from the 
saddle-point in the direction of the stripe. This soften- 
ing is absent at lower doping [cf. Fig. [TO](n/! = 0.0625, 
rif, = 0.083)] and indicates the susceptibility towards a 
magnetic instability along the stripe at higher doping. 
It has been argued L73i] that this roton like minima ex- 
plains INS in untwinned samples of YBCO by Hinkov 
and co-workers i6lll . 

Finally, at higher energies the excitations again 
spread out and form a ring-shaped feature around Qaf 
with a small anisotropic intensity modulation (Fig. 
|8] upper panel). Note that our spectra are for one- 
dimensional stripes without any orientational average so 
that C4 symmetry is broken. The resulting quasi two- 
dimensional excitation spectra above the saddle-point 
energy are in marked contrast to the one-dimensional 
characteristics obtained from localized spin models 1741 
7511 ■ This has profound implications on the interpreta- 
tion of INS data from detwinned samples j 68 , 6^ which 
also show pronounced 1-D characteristics at low but a 
2-D magnetic distribution of magnetic energies at high 
energies in accord with our result. 

The high energy distribution of magnetic intensity 
can be understood from the energy weighted sum rule 



Jo 



dajojx'^iaj) = -— 2j(ek-ek+q)nk,<r(5) 



where ek denotes the single-particle dispersion of the 
underlying itinerant model and nk.u is the number op- 
erator for particles with momentum k and spin cr. In 
the presence of time reversal symmetry and inversion 
("k.o- = n-k-cr) this can be simplified to 



(6) 



and 6 runs over the vectors connecting a specific lattice 
site to its respective neighbors with only one represen- 
tative of the pair (S, -6). Thus in case of a square lattice 
6 - x,y for nearest neighbors, 6 = x-i-y, x-y for 
next nearest neighbors etc. and Tg denotes the kinetic 
energy in J-direction. Since weights more the exci- 
tations at higher energies it turns out from Eq. (|6]l that 
with increasing oj the intensity in the transverse mag- 
netic susceptiblity should be confined to wave-vectors 
around Qaf - 

Consider now the case of stripes oriented along the y- 
direction. Although we are dealing with a perfectly or- 
dered quasi ID structure we find from our calculations 
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that e.g. for d — 4 BC stripes the kinetic energy along 
the y-direction is only ~ 10% larger than that perpen- 
dicular to the stripes. Neglecting the contribution from 
next-nearest neighbor terms Eq. ^ becomes 

M'^^sin^f^l + ^sin^f^) 

so that even in the presence of stripes the high energy 
magnetic response is essentially two-dimensional but 
slightly anisotropic around Qaf- In Fig. [8]this is exactly 
the distribution of the optical magnetic excitations one 
finds above the saddle-point energy. We anticipate that 
this feature can also be observed in the spin glass-phase 
of LSCO where due to the imbalance of twin domains 
the one-dimensional character of low energy spin exci- 
tations has been resolved by inelastic neutron scattering 
experiments (6^. Above the measurements seem 
to indicate that the excitations acquire again the (two- 
dimensional) character of the magnons in the undoped 
system in agreement with the prediction for a striped 
ground state [29] . 

Another interesting aspect of the sum rule Eq. |5] that 
has passed unnoticed so far, is that it allows to measure 
the change in kinetic energy on entering the supercon- 
ducting state in a way that is alternative to the proposal 
in Ref. |63]. The present proposal needs accurate mea- 
surements in the higher part of the spectra which should 
become available with the new pulsed neutron sources. 

The overall doping evolution of the vertical stripe 
magnetic excitations together with availabe experimen- 
tal data is depicted in Fig. [lOl Here we show the exci- 



tations along scans (a,b) sketched in Fig. |9] It should 
be noted that the magnetic excitations for SC stripes are 
almost identical with only slight differences regarding 
the intensity distribution and the gap structure of opti- 
cal branches i72ll . It turns out from Fig. [10] that for 
doping «/, > 0.125 the saddle-point excitation at 
Qaf rapidly shifts to higher energies together with an 
increasing softening of the magnetic excitations along 
the stripe. In fact, the value of a>s is determined by 
the magnetic coupling across the domain wall. Because 
this coupling is mediated by virtual hopping processes 
of holes between stripe sites and adjacent AF regions 
it strongly depends on the stripe filling factor v as de- 
fined above. As we have shown above for «/, < 1/8 the 
filling of the core of the stripe remains practically con- 
stant and cl>s has a weak dependence on doping. On the 
other hand for «/, > 1 /8 doping proceeds by changing 
the filling of the stripes at constant incommensurability. 
This produces a rapid renormalization of the effective 
exchange interaction across the core and a concomitant 
rapid increase in ojs- 

In Fig. [TT]we compare the doping dependence of oJs 
with available experimental data |56, 57, 59, 60, 62, 64ll . 
The latter seem to indicate lf59ll that in LSCO the saddle- 
point excitation at Qaf stays at least constant beyond 
«/, - 1 /8 in contrast to what is obtained from the spectra 
of doped stripes shown in Fig. (TO] It is conceivable that 
instead of doping additional holes into the domain walls 
these (fluctuating) textures in LSCO stay half-filled be- 
yond «/, = 1/8 leading to phase separation between 
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Figure 1 1 : Doping dependence of the saddle-point excitation tOj at 
Qaf for diagonal site- and bond-centered stripes and vertical bond- 
centered stripes. Circles with error bars are experimental data from 
Refs. 57. 59, 60, 62, 64]. The horizontal dashed line starting from 
= 0.125 is the expected behavior for a phase separation scenario 
for stripes. 



free carriers and those which are involved in the stripe 
correlations |25]. On the other hand Wakimoto and 
collaborators 1 24] find that at large doping «/, - 0.3 scat- 
tering starts around 80meV which is compatible with a 
strong increase of u)s as we predict for the doped stripes. 
Thus more experimental work is needed at intermediate 
dopings to clarify the situation. 

With underdoping the computations increasingly 
overestimate the energy of the experimentally deter- 
mined cOs, especially in the diagonal phase. The main 
reason for the disagreement in this doping range is prob- 
ably the increasing disorder character of the stripes. To 
support this idea we have performed 1I29I1 linear spin 
wave theory calculation of disordered diagonal bond 
centered stripes. Is is found that already a small amount 
of disorder leads to a softening and broadening of 
which would bring the theoretical data in Fig. [TT] closer 
to the experimental ones. 

Recent experimental effort has been made to measure 
magnetic excitat ions in LSCO with resonant inelastic 
x-ray scattering f 164"]. Because this technique essen- 
tially measures in the nuclear Brillouin zone (i.e. close 
to momentum q = 0) we show in Fig. [12] the corre- 
sponding spectra for d - A y-axis oriented bond cen- 
tered stripes. Consider first the scan starting from mo- 
mentum q - (0, 0) [labeled (d) in Fig. |9l in the direction 
along the stripe. Because the system is still AF ordered 
along the stripe (cf. Fig. |2| one has a doubling of the 
unit cell in this direction. At the same time the doping 
of the stripe is approximately v » 0.5 so that the low 
energy structure of the spectra reflects the continuum of 
spinon Uke spin-flip particle-hole excitations for the ac- 
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Figure 12: Magnetic excitations in the nuclear Brillouin zone for d = 
4 y-oriented stripes. The spectra correspond to the scans labeled (c,d) 
inFig.|9] 



tive band shown in Fig. |6] At higher energy, part of 
the (high intensity) optical magnon band overlaps with 
the spin-flip continuum and dominates the spectrum. In 
the direction perpendicular to the stripe scan (c) is con- 
nected to scan (a) by the addition of a reciprocal lat- 
tice vector (37r/4,7r). Therefore the low energy magnon 
band which disperses from (0, 0) to (n, 0) is the replica 
of the Goldstone mode which has been discussed above. 
The same holds for the higher energy branches which, 
however, in the nuclear and magnetic zones strongly dif- 
fer in weight. 



7. Optical conductivity 

In the previous sections our considerations where 
mostly based on the extended one-band model, which 
should be appropriate as long as interband transitions 
(as in the case of spin excitations) do not play a sig- 
nificant role. In the present section we discuss the 
optical conductivity of striped systems which involves 
also higher energy excitations. We therefore investi- 
gate the corresponding spectra in the framework of the 
three-band model which is again treated within the time- 
dependent Gutz wilier approach ]'25','2^. 

It is also possible to obtain an optical conductivity 
that in the main features resembles experiment, using 
a one-band Hubbard model. Indeed, one could cer- 
tainly argue from the inset to Fig. [TJ) that the transi- 
tion between LHB and UHB mimics the charge-transfer 
excitation of cuprates and therefore should provide a 
reasonable starting point. However, the point is that 
the charge excitations will also be strongly influenced 
by the respective stabilitity of bond- and site-centered 
stripes which are energetically degenerate in the one- 
band model. This leads to soft lateral fluctuations of the 
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stripes (phason mode) which are optically active. Be- 
cause of the degeneracy between the two textures, the 
phason is soft in the one-band model at low doping. 
In contrast, in the more realistic three-band model the 
two textures become only degenerate at higher doping 
which, as shown below, is in good agreement with ex- 
periment. 




Figure 13: Optical conductivity of the three-band Hubbard model 
evaluated within the time-dependent GA for various hole concentra- 
tions. In the inset: The evolution of the mid-infrared peak as a func- 
tion of doping. The Drude contiibution has set to zero in this case. 

Based on the three-band hamiltonian Fig. [T3] reports 
the optical conductivity for AF and inhomogeneous 
ground states at various dopings. The parameter set has 
been taken from the first principle computations of Ref. 
1 1671] so we have no adjustable parameters. 

At very dilute doping («/, < 0.03) due to the long- 
range Coulomb interaction (not included in our calcula- 
tions) each hole will be close to an acceptor preventing 
the formation of stripes. The RPA optical conductivity 
for the corresponding single-hole solution («/, = 0.028) 
leads to the formation of a doping induced MIR band 
close to 0.5eV and doping induced transfer of spec- 
tral weight from the charge transfer (CT) band to the 
MIR regi on in agreement with experiment in this dop- 
ing range OlSlll . 

For distant stripes (d = 7) the single-hole MIR band 
now splits into two bands. The one at higher energy is 
a band of incoherent particle-hole excitations close to 
1.3eV which provides a theoretical explanation for the 
shoulder to the CT excitation seen at the same energy in 
optical absorption through LSCO thin films 1 165) , and 
electron energy loss spectroscopy lll66ll . 

The low energy MIR peak is a collective mode asso- 
ciated with the soft lateral displacements of the stripe 
mentioned above. It is located at 0.3 eV at low do- 



ing and becomes soft only at optimum doping where it 
merges with the Drude response in good agreement with 
experiment. This soft collective modes are low energy 
excitations which have no counterpart in a noninteract- 
ing system breaking the paradigm of Fermi liquid the- 
ory Ii26il . It is also possible that coupling to these modes 
contributes significantly to the pairing. It is worth to 
remark that new developments in time resolved spec- 
troscopy show [169,1 that the high energy excitations be- 
low the CT band do not participate in the pairing but 
excitations at the CT band do. The region below 1.6eV 
has not been measured yet. 

It is interesting to remark that also in nickelates the 
MIR peak has been attributed to the for mati on of mid- 
gap states due to the stripe ground state lll70ll . 

8. Dynamical stripes 

The investigation of stripe 'dynamical properties' in 
the previous sections concerned the excitations on top 
of a stripe ground state with long range order. However, 
as discussed in Secs. l2l6l in most cuprate materials long 
range static spin correlations are not observed and also 
direct evidence for static charge order is quite elusive. 

Local probes like STM often point to a broken trans- 
lational symmetry in real space specially when the high 
energy part of the spectra is observed. Momentum space 
probes like ARPES and INS resemble those of an or- 
dered state at high energy but interpolate to a state with- 
out broken symmetry at low energy. This points towards 
a dynamical or a glassy nature of the stripe correlations. 

There are two ways in which this dichotomy be- 
tween low and high energy quasiparticle excitations can 
be solved: 1) proximity to a quantum/classical critical 
point in the disordered side (Sec. 18. 11 1. 2) Stripes with 
long range order but with protected low energy quasi- 
particles (Sec. l8.2l i. For the last possibihty there are also 
two variants: 2. a) Long range order in spin and charge 
and 2.b) Long range order only in the charge. 

8.1. Proximity to a quantum/classical critical point 

The classical version of scenario 1) consists of ther- 
mally disordered stripes at finite temperature. This 
has been the route followed by Vojta and collaborators 
1 781] who have investigated a Landau model for coupled 
charge and spin fluctuations in the Born-Oppenheimer 
approximation. Although not mentioned explicitly, this 
computation also describes a glassy stripe state. 

The quantum version of scenario 1) follows from this 
qualitative argument: In physical dimensions a system 
may have long (but finite) ranged order parameter spa- 
cial correlations which are also long lived close to a 
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quantum critical point. This defines a fluctuating fre- 
quency cl>o above which the systems appears to be or- 
dered (or at least critical). Then, for energies larger than 
cjQ with respect to the Fermi level, the spectral function 
should resemble the spectral function of an ordered sys- 
tem. This high energy spectral weight which normally 
one would term "incoherent" may in reality carry im- 
portant information on the momentum structure of the 
close-by ordered phase. This momentum structure may 
be determined with ARPES from the momentum dis- 
tribution Hk by integrating the spectral function over a 
broad energy window or by STM and INS experiments 
by direct examination of the high energy spectra. 

On the other hand, electrons at lower energies aver- 
age over the order parameter fluctuations and "sense" a 
disordered system. In this limit we expect Fermi liquid 
quasiparticles with all their well known characteristics 
like a Luttinger Fermi surface. This dichotomy in the 
momentum structure of low and hig h energy q uasipar- 
ticles is consistent with ARPES data jlOdllOTll on lan- 
thanum cuprates. 

In a series of papers 11171 , 



a) 



172, 173 



17411 which 

were based on a Kampf-Schrieffer type approach f 176ll 
we have worked out the scenario of fermionic quasi- 
particles coupled to dynamical charge order fluctua- 
tions with regard to several experimental observations. 
These include the ang ular dependence of the quasiparti- 
cle weight in Bi2201 lll77ll . the isotope shift of the high- 
energy electronic dispersion iflTStl . and t he dichoto my 
in the Fermi surface of high-Tc cuprates 1 106l 107 ] as 
mentioned before. 

The Kampf-Schrieffer approach describes quite suc- 
cessfully situations in which the quasiparticles are cou- 
pled to locally well formed order-parameter fluctuations 
but can not be derived from a microscopic approach. On 
the other hand a formal derivation close to a quantum 
critical point leads to quasiparticles coupled to diffusive 
collective modes (CM)s. This leads to a form which is 
customary in quantum critical p heno mena and has often 
been used for spin fluctuations il79ll 



1 



- (8) 



where A - c,s refers to charge or spin CMs substan- 
tially peaked at characteristic wavevectors ^ (tt, n) 
and qc ~ (+;r/2, 0), (0, +7r/2) respectively. Here the 
mass m,i is the minimum energy required to excite the 
CM and va is a fermion scale setting the CM momen- 
tum dispersion. This dispersion is limited by an energy 
cutoff A^. The dimensionless quantities m^/A^ are the 
inverse square correlation lengths (in units of the lattice 




Figure 14: Constant frequency cuts of the electronic dispersion for 
dynamical d = 4 bond centered Gutzwiller stripes, (a) a> = Ef, (b) 
a) = -0.6f. Parameters: U/t = 8, f /t = -0.2. The frequency structure 
is obtained from Eq. ll2l with ojo = O.lt and v- = 0.02. 



spacing), which measure the typical size of ordered do- 
mains. The ica term establishes the low-energy diffusive 
character of these fluctuations due to decay into particle- 
hole pairs, whereas above the scale set by Q.^ the CM 
has a more propagating character This approach has 
been used to investigate the crit ical b ehavior of the low- 
frequency optical conductivity f ISO'] and, more recently, 
to calculate the Raman respons e fun ction including self- 
energy and vertex corrections lllSlll . In particular, as- 
suming as mediators spin and charge fluctuations with 
different characteristic wavevectors, it was possible to 
selectively individuate the charge and spin contribu- 
tions in the Raman response function. By fitting the 
experimental spectra in La2-xSrxCu04 single crystals 
within the above theoretical framework, it was found 
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that both charge and spin fluctuations were contribut- 
ing to the quasiparticle scattering (a distinctive feature 
of stripe fluctuations). However, it was also found that, 
upon increassing doping, charge CMs progressively ac- 
quire more relevance while spin fluctuations become 
weaker and weaker This indicates that upon increas- 
ing doping into the overdoped regime, stripes naturally 
evolve into harmonic damped charge-density-wave fluc- 
tuations, while spin modes eventually loose weight. 

8.2. Protected low energy quasiparticles 

STM il3d 131 1 experiments provide strong evi- 



dence for the simultaneous existence of low energy (Bo- 
goljubov) quasiparticles and manifestations of truly bro- 
ken translational symmetry at h igh e nergy. 

In another recent approach 111841 174] we have de- 
veloped a phenomenological description with broken 
translational symmetry but a dynamical order parame- 
ter which applies to this situation. The theory also ac- 
counts for the contrast reversal in the S TM spectra be- 
tween positive and negative bias 117511, in contrast to 
mean field like static order where this effect in general 
occurs away from the Fermi level. Here we show that 
this theory can be combined with the GA computations 
discussed in Sec. |4]in order to phenomenologically ac- 
count for the electronic structure of a fluctuating stripe 
state. 

The GA is a renormaUzed mean-field theory which 
allows for the definition of an effective Gutzwiller 
hamiltonian (cf. e.g. Ref. Li34J) 



ijo- io- 



(9) 



where the c^" are now (creation) annihilation operator 
of the Gutzwiller quasiparticles. The are renormal- 
ized hopping amplitudes and /l,o- correspond to local 
chemical potentials. Both parameter sets depend on the 
value of the onsite repulsion but also on the charge and 
spin structure of the inhomogeneous ground state. 

We want to construct an ansatz for the self energy due 
to electron-electron interactions. In order to ensure that 
the self-energy is physical (i.e. that obeys all analyti- 
cal properties of an electronic self energy) we map the 
problem to that of electrons coupled to a set of auxil- 
iary states. The idea is to enlarge the Hilbert space by 
introducing the coupling to an additional set of local- 
ized states (described by creation and destruction oper- 
ators flj^ which thus results in the following (auxiliary) 
Fano- Anderson Hamiltonian, 



As a result the Greens function for the Gutzwiller 
quasiparticles depends on an additional self-energy 
which is determined by the electronic structure of the 
/-states via 



2,(w) = y + A,- = Ucj) + A/. 



(10) 



with the constant A, controlling the high frequency 
limit. 

Note that this kind of approach automatically main- 
tains the correct analytical properties of the (frequency 
dependent) self-energy which replaces the static single- 
particle potential in Eq. |9] 

For the frequency structure in Eq. \W\ we choose a 
two-pole Ansatz 



fi(oj) = 
A/ = 



Pi 



(jj - Li>0 10 + COq 



(11) 
(12) 



with Q', H- /?, - and controls the strength of the 
scattering. Thus at each site an asymmetry is intro- 
duced in the spectral distribution which below is taken 
to be proportional to the charge modulation. The con- 
stant A; guarantees the limit S,(w = 0) = at each 
site. Denoting by A the average local chemical po- 
tential the weights of the poles are choose as a, — 
1/211 -tanh(l -^,/I)]. 

Fig. [14] reports the corresponding scans of the elec- 
tronic structure at the Fermi energy (a) and at an en- 
ergy of 0.6f below Ef. Clearly the low energy electronic 
structure resembles that of a homogeneous system due 
to die fact that the Ansatz Eq. [T2]leads to E,(w = 0) = 0. 
Instead at high energies w » wo the scattering is con- 
trolled by Z,((y — > oo) = A, which according to the 
above definitions is governed by the local charge and 
spin structure. One can clearly see the similarities be- 
tween the electronic structure ai o) - -0.6f (cf. Fig. 
[T4b ) and the reconstructed Fermi surface from the static 
stripe solution shown in Fig. |7^. 

Therefore this phenomenological theory reproduces 
our scenario of low energy protected quasiparticles but 
emerging spectral properties of charge and spin order 
at large e nergi es. In contrast to the Kampf-Schrieffer 
approach 117611 which is based on homogeneous ground 
states the present theory explicitely describes systems 
with broken symmetry, however, this broken symmetry 
manifests only at high energies. 



16 



9. Conclusions 

In this paper we have reviewed evidence for 
stripe correlations in cuprate superconductors based on 
Gutzwiller variational calculations supplemented with 
Gaussian fluctuations. Because our aim is to pro- 
vide also a quantitative analyis of static and dynamical 
properties, the appropriate parameter set for lanthanum 
cuprates based on the extended Hubbard model was de- 
rived in Sec. [3] in some detail. O ur val ue of U/t - 8 
agrees with estimates from Ref. 11 14411 and the next- 
nearest neighbor hopping t' /t - -0.2 is in the ra nge o f 
values as derived from LDA computations in Ref. 18511 . 

The main emphasis of this review is put on the spin 
excitations on top of stripe ground states. As we have 
discussed in Sec. |6]the time-dependent Gutzwiller ap- 
proach yields an excellent quantitative agreement with 
the spectra observed by INS when the stripes are static, 
i.e. in the LBCO compounds. Upon underdoping and 
in the non-codoped system our calculations in general 
overestimate the excitation energies. We attribute this 
discrepancy to the much reduced correlation length as 
extracted from the low energy spectra due to the dynam- 
ics of the stripes and disorder eff'ects [29J. In this regard 
it would also be interesting to compute the spin exci- 
tations on top of the dynamical stripe solutions which 
have been discussed in Sec. |8] 

As one approaches the overdoped regime the nature 
of stripe correlations in the ground state is not so clear 
since one has to take into account the melting of the 
stripes and non-Gaussian fluctuation eff'ects will be- 
come important. Also it is conceivable that phase sepa- 
ration between the underdoped stripes and an overdoped 
Fermi Uquid occurs as discussed in Ref. ll25ll . Evidence 
for this scenario was presented in Sec. |6]from the dop- 
ing dependence of the saddle point energy. Due to the 
fact that in the phase separation scenario the population 
of the domain walls by holes saturates a further increase 
of the saddle point energy with doping is harnpered in 
agreement with the INS experiments in Ref. 115 911 . On 
the other hand, Ref. fl^ is in accord with a substantial 
increase of the saddle point energy. 

Another possiblity is that the anharmonic structure 
of the stripes is gradually lost and stripes become har- 
monic dynamical and damped charge-density wave fluc- 
tuations. This leads to the more standard scenario of a 
second order transition ending into a quantum critical 
point slightly above optimal doping[3] driven by a frus- 
trated phase separation mechanism[4|]. In this case the 
combined effect of charge and spin fluctuations can sim- 
ply be treated within perturbation theoryl 186|. 

The natural question which arises is whether the re- 



sults of Sec. |6] may then also be relevant for an un- 
derstanding of spin excitations in YBCO. In this com- 
pounds the acoustic dispersion shows a low-energy spin 
gap (A ~ 30meV at optimal doping) and thus does not 
reach the incommensurate wave-vector at w = 0. There- 
fore it is more likely that the system is in a quantum 
disordered s pin phase as suggested by the ladder theo- 
ries 1 74, 751176^. Alternatively a scenario of fluctuating 
stripes where also the charge loses its long-range order 
can capture the effect of a spin- gap . ifTill One expects 
that the systems show short range order with a corre- 
lation length of the order of A/ (he) ~ 5a and that for 
energies larger than A the system resembles that of an 
ordered phase. For this reason one can expect that a 
computation as the one discussed in Sec. |6]is suitable 
for a description of the universal high e nergy spin re- 
sponse (sE, '58', '59', S lil P, '6?, '66 67', '69^ . Again, to 
capture the excitations over the full energy range would 
require an approach which incorporates the dynamical 
nature of stripes as the one proposed in Sec. [8] 

Concluding, we have shown that the analyis of nu- 
merous experiments provides strong evidence for stripe 
correlations in high-T^ superconductors. This more 
or less hidden electronic order can account for both 
normal-state anomalies and high-temperature supercon- 
ductivity in the cuprates. In fact, the proximity to an 
instability (particularly if it is a second-order "critical 
line" ending at zero temperature into a quantum critical 
point) marking the onset of order naturally brings along 
abundant fluctuations and leads to strongly temperature- 
and doping-dependent features, which account for the 
non-Fermi liquid properties and for a strong pairing in- 
teraction. The participation of such stripe fluctuations 
in the 'pairing glue' is sugge sted by the recent analy- 
sis of Raman scattering data fl86'l but in any case ad- 
ditional work has to be done in order to elucidate the 
source of strong scattering/pairing between the charge 
carriers leading to high-T^. 
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